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ABSTRACT: 

In this paper, flow of air through a pressure regulating and shut-Off Valve (PRSOV) has been 
solved numerically with an objective to determine the pressure of the PRSOV pilot chamber. Flow 
through the valve is transient, compressible and turbulent in nature. Flow has been solved using 
ANSYS FLUENT coupled with a special user defined function (UDF). UDF is used for dynamic 
meshing, automatic data acquisition from ANSYS FLUENT solver and to schedule the numerical test 
procedure. Within the UDF, valve inlet pressure is varied in a stepwise manner. For every value of 
inlet pressure, transient analysis leads to a quasi-static flow through the valve. Spool forces are 
calculated based on different pressures at inlet. From this information of spool forces, pressure 
characteristic of the pilot chamber of passive control circuit has been derived. The same 
characteristics have also been obtained after modelling the flow analytically. Both the results have 
been compared. It is observed, CFD analysis of the flow has led to improved results. 

KEYWORDS: Pressure regulating and Shut-Off Valve, pilot flow, dynamic meshing, flow force, 
pressure regulation, Force Balance, CFD 



I. INTRODUCTION 

Pressure regulating and Shut-off valve is a critical component used in many fluid power systems, 
especially in aircraft. The valve is used in Environmental Control System of the aircraft. The valve is installed in 
a pipe line in an air flow path where available inlet pressure to the valve may vary within a wide range. The 
main function of the valve is to regulate the air flow and maintain almost constant pressure at its outlet. 
Schematic diagram of the valve is shown in Fig. 1. 

A passive control circuit (Pilot chamber pressure control circuit in Fig. 1) is required to achieve the 
above objective. Depending upon the valve inlet pressure, this passive control circuit sets the criterion for 
position of spool inside the valve. Determination of characteristics of the control circuit is a crucial part in 
designing the control circuit and hence the pressure regulated valve. This requires complete information of flow 
forces acting on the valve spool. Calculations based on theoretical model can provide some initial approximate 
information. With the advent of numerical power, advanced, sophisticated and reliable CFD tool is very useful 
in determining the various flow forces acting on the spool. 

Leutwyler and Dalton [1] showed the potential of CFD tool in analyzing compressible, turbulent flow 
through butterfly valves. ANSYS FLUENT was used for the analysis. Amirante et al. [2-3] evaluated flow 
forces on an open centre direction control valve. Amirante et al. [4] again modelled direct proportional valve 
using CFD. Compensation techniques based on spool profiling were used to balance the flow force at different 
level of valve openings. Chen et al. [5] reported flow visualization using CFD in a ball valve. Chattopadhyay et 
al. [6] have investigated turbulent flow structure inside a Pressure Regulating and Shut-Off Valve using ANSYS 
FLUENT. In this work, both 2-D and 3-D simulations were performed with a conclusion that 2-D model could 
predict the flow coefficients satisfactorily. Song et al. [7] have reported 2D dynamic simulation of a pressure 
relief valve using CFD. 
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Orifice Solenoid valve 




VIEW A 

Fig. 1. Schematic diagram of pressure regulated valve 
Variable Orifice opening — ► — 



Pilot Chamber 




View A 

In this study, transient, compressible and turbulent flow through a pressure regulated valve has been 
solved using commercial software ANSYS FLUENT [10]. As the valve inlet pressure changes, spool inside the 
valve as shown in Fig. 1 must move in appropriate direction to keep outlet pressure constant. Starting from the 
instant of inlet pressure change, a transient analysis is performed to reach the state of quasi-static flow through 
the valve. 

II. ANALYTICAL METHOD 

Pressure Regulating and Shut-Off Valve (PRSOV) consists of valve body, spool and pilot chamber. 
Inlet pressure (Pi) of the valve may vary from 207 kPa to 758 kpa. Main function of the valve is to keep the 
outlet pressure (P 2 ) at almost constant value. Spool inside the valve body is subjected to a number of forces. 
When there is a change in inlet pressure, spool moves in appropriate direction and outlet pressure should remain 
almost constant value. 

2.1 Force balance on the spool 

spool moves under the action of four forces, F b F 2 , F 3 and F 4 as shown in Figure 1. Force F x arises due 
to the pilot chamber air pressure (P 0 ). Pilot chamber air pressure is controlled by datum setting of relief valve as 

shown in Fig. 1. The effective area on which this pilot pressure acts is 7id 0 2 /4. Therefore, ¥ 1 = P 0 x — and acts 

from left to right direction as shown in Figure 1. 

As the high pressure air flows through the variable orifice, air pressure reduces from inlet pressure (Pi) 
to outlet pressure (P 2 ). This outlet pressure acts on the other side of the spool head. Thus a force (F 2 ) acts on the 

spool from right to left direction as shown in Fig. 1. Here, F 2 = P 2 X 
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Air pressure changes from Pi to P 2 as it flows through the orifice. Orifice is formed by spool flange and 
the valve body thickness (t 0 ). Since t 0 is small enough, so pressure variation is assumed linear. Therefore, 
average pressure (Pi +P2)/2 acts on spool flange over an annular area with diameter (d 0 + t 0 ) and thickness t G . 

Here the force, F 3 = ^ 1+ ^ 2 ^ x Ti(d 0 + t 0 ) X t 0 . Since t 0 is small value, so higher order term of t 0 can be 

neglected. Then, F 3 = ^ 1+ ^ 2 ^ x Tid 0 X t 0 and the force (F 4 ), which acts on similar area on the other side of the 
spool flange, is given as F 4 = P x X Tid 0 X t 0 . 

Spool remains at equilibrium state under the action of these four forces. Therefore, ¥ 1 + F 3 = F 2 4- F 4 

Or,P o X^ + ^XTTd 0 Xt 0 =P lX 7ld 0 Xt 0 + P 2 X^ 

or, P 0 = X P x + fl - X P 2 . 

0 2d 0 1 V 2d 0 J Z 

III. CFD METHOD 

3.1 Mathematical modelling 

The ANSYS FLUENT package has been used to solve the axisymmetric flow in the domain shown in 
Fig. 2. The partial differential equations for the conservation of mass, momentum, energy, turbulence kinetic 

energy ^ and turbulence dissipation rate 8 along with the algebraic equation for the perfect gas laws have been 
solved. As given below, these involve time t as the independent variable, the dependent variables of density P , 

velocity v , pressure p , temperature T and the flow parameters of kinematic viscosity V 1 , turbulent viscosity 



A/r dp/dt + V .(pv) = 0 

Mass conservation equation: r xr 

lation: 

d(p7)/dt + V .(pvV) = -V j9 + p(v / +v f )V.{Vv + (Vv) r -(2//3)V.v] 



1 , molecular diffusivity 1 and turbulent diffusivity 1 . 

[ass conservation equation: ^P^ [ 
Momentum conservation equation: 

d(pv)ldt + V .(/?vv) = -V j9 + 

Energy conservation equation: 

dT ldt + V.VT = -(a l + a t )V 2 T + (VV):[p(v l +v t ){VV+{VV) T -(21 
Turbulence Kinetic Energy Equation: 

d(pk)ldt + V (pvk) = V .[p(v l + v t l<J k )Vk]+G k - ps - 2pskl(yRT ) 
Turbulence Dissipation Equation: 

d(ps )/dt + V .(pvs) = V .[p(v l +v t /cr k )V s]+ pC x Ss- pC 2 s 2 l(k + 4vs ) 
where ^=[{ Vv + (Vv) r }:{Vv + (Vv) r }/2] 1/2 ^ 

G k=v t s 2 

C 1 and C 2 are constants and ° k and ° 8 are turbulent Prandtl number respectively, R is the gas constant and 

^ is the specific heat ratio. Consistent with the realizable k-s model due to Shih et al. [13], the equations and 
the values of the different parameters are chosen as 

v t =C k 2 Is 



a t= v r v t / ( a i Pr f) 
and 



ju=ju o (T/T o ) 3/2 {(T o + s)/(T + s)} 



where ^ 0 and T ° are reference viscosity and temperature, s is Sutherland constant and Pr 1 is the turbulent 
Prandtl number, which have been taken respectively as equal to 1.716 x 10-5 Pa-s, 273.11 K, 110.56 K and 0.85. 
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3.2 Boundary Conditions 

Following boundary conditions were used for this analysis. 

1. Inlet boundary where a fixed pressure condition is used. The inlet temperature is fixed at 633 K. 

2. As discussed above, outlet boundary of the flow domain is open to ambient. So, pressure and 
temperature have been set to 101325 Pa and 298.15 K, respectively. 

3. No-slip boundary condition is assigned for all the walls. The valve wall was assumed to be at adiabatic 
condition. 

4. For prescribing turbulent quantities at the boundaries, several options are available. We have used a 
condition prescribing the level of turbulent intensity (TI) at around 5-10% in the incoming fluid stream 
which is a reasonable level of value used by researchers. Hydraulic diameter of 47.5mm in the 
incoming fluid stream has been considered. 

The initial conditions are the steady state result corresponding to the pressure p x at the inlet of the valve 
body equal to 207 kPa and all other boundary conditions as stated above. 

Flow geometry has been approximated with a 2D axisymmetric domain as shown in Fig. 2. 
Unstructured quadrilateral mesh has been generated using commercial software package ANSYS ICEM CFD. 
Dynamic mesh motion along with the movement of spool boundaries has been considered during mesh 
generation. Entire fluid flow domain has been divided into two parts as shown in Fig. 3. Mesh for zone 1 does 
not possess any motion, whereas mesh for zone 2 undergoes dynamic motion. Wall boundaries comprising spool 
surfaces move depending upon the conditions as mentioned in FLOWCHART in Fig. 5. Some of the internal 
edges also move with spool surfaces. One special internal edge has been selected to move with spool surfaces 
but at exactly half the speed of spool movement. This special movement serves a very important purpose; finer 
boundary layer of the variable orifice remains intact though the spool opening reduces. Initial mesh contains 0.2 
million quadrilateral cells approximately. 

v f^ ' 1 

Fig. 2. 2D axi- symmetric flow domain 




Fig. 3. Flow domain showing two different zones of fluid 
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Fig. 4. Structured quadrilateral mesh around orifice opening 



3.4 Methodology of calculation 

Valve inlet pressure may vary from 207 kPa to 758 kPa, whereas valve outlet pressure remains within 
200 kPa to 241 kPa. Air flow further downstream of valve outlet is approximated with the help of an orifice 
flow finally open to ground level condition. Calculation domain, therefore, includes the valve geometry along 
with the appropriate pipe length before and after the orifice. 2D axisymmetric flow geometry is shown in Fig. 2. 

Transient flow through this flow domain is solved in ANSYS FLUENT to determine various spool 
openings depending on various inlet pressures. Now, when the inlet pressure is at the minimum value, i.e., 207 
kPa, spool opening should be 100%. Therefore, solution of this flow sets the initial condition for transient 
calculation. In transient calculation, inlet pressure is varied in a step wise manner. This is done via a User 
Define Function (UDF). In reality, when inlet pressure is changed by a small amount, flow through the valve 
will take some time to die down the transient phenomena. In CFD, this quasi-static state is achieved after 
solving the transient flow for a number of time steps. For each time step, the flow convergence is checked using 
UDF that allows the calculation to move to next time step only after meeting the required criteria for 
convergence. Also, a check, whether quasi-static condition has been achieved, is carried out with the help of 
same UDF before proceeding to next time step. This UDF is described with the help of a FLOWCHART given 
in Fig. 5. 

When the transient dies down pressure at valve outlet is checked for the desired outlet pressure. If 
desired pressure is not obtained, then spool opening is reduced by a delta amount and further transient 
calculation is carried on. This whole process is automated using the UDF. When the desired pressure is 
obtained, UDF directs ANSYS FLUENT to write the Case and Data file for further post processing. Spool 
opening and valve outlet pressure for each time step are written in a DAT file with the instruction written in 
UDF. This file is required for further reference. After achieving the desired pressure at valve outlet, UDF directs 
ANSYS FLUENT to increase the valve inlet pressure by a small amount and again the whole loop of operations 
repeated. 
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Fig. 5. FLOWCHART for special UDF 

IV. RESULTS AND DISCUSSION 

Flow through the pressure regulated valve is governed by variable orifice flow. Fig. 6 shows the spool 
opening in mm against the gauge pressure at valve inlet. It is observed here that spool opening reduces by 86% 
as valve inlet pressure increases from 240kPa to 360kPa. Thereafter, the gradient of change in spool opening 
with inlet pressure is very slow. This happens because pressure is almost fixed at valve outlet. 

Fig. 7 shows the variation of mass flow rate with the change in inlet pressure. There is no significant 
change in mass flow rate. This is justified as two things happen simultaneously. 




240 280 320 360 400 440 480 520 560 600 640 680 720 760 
Inlet Pit ssii ir i UP. 1 1 

Fig. 6. Spool opening plotted against valve inlet pressure. 
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Fig. 7. Variation of mass flow rate through valve with change in inlet pressure. 
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Resistance in flow path increases as spool opening reduces, which would try to decrease mass flow 
rate. However, the increase in inlet pressure almost nullifies the above effect. 

In Fig. 8 flow forces acting on spool surfaces are plotted against the inlet pressure. The force is acting 
axially from right to left direction with respect to Fig. 1. Both analytical and CFD results are plotted here. In 
both the cases, it is observed that the net force increases with increase in the inlet pressure. This owes to the fact 
that the spool opening creates an orifice and the air pressure acting on the orifice side of spool flange is less than 
the pressure that acting on the outer side of the spool flange. The derivation in Section 2.1 makes this aspect 
clearer. Higher absolute force is predicted by the CFD analysis. This pressure distribution and the consequent 
prediction of spool force is of course more realistic in view of the more detailed nature of the CFD analysis in 
comparison to a lumped calculation. 

As the spool is subjected to a flow force from right to left, a counter force is required to keep the spool 
at the desired location. This force comes from pilot chamber air pressure. Air pressure in the pilot chamber is 
maintained by a Pilot Chamber Pressure Control circuit. There is a spring loaded Pressure Relief Valve in that 
circuit. The required pilot pressure to balance the spool is calculated and plotted in Fig. 9. Both Fig. 8 and Fig. 9 
show that CFD analysis has led to sufficient improvement of the pilot chamber pressure characteristic with 
respect to the analytical results. 
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Fig. 8. Variable flow forces on spool surfaces with inlet pressure variation. 
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240 280 320 360 400 440 480 520 560 600 640 680 720 760 

Inlet Pressure (kPa) 
Fig. 9. Pressure characteristic of passive control circuit. 



Fig. 10 shows the contour plots of static pressure at 100%, 60%, 30% and 10% spool opening. 
Corresponding inlet pressures are 241 kPa, 255 kPa, 276 kPa and 414 kPa respectively. A low pressure zone is 
created at immediate downstream of variable orifice. After that pressure recovery takes place. The size of low 
pressure zone increases with the increase in valve inlet pressure. This becomes clear from the stream function 
plots in Fig. 11. The plots clearly reveal the recirculation bubble and the associated downstream flow divergence 
in each case. Due to the increase in the size of the recirculation bubble with increase in the inlet pressure, the 
pressure recovery zone is pushed further right. 
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30% opening 10%opening 

Fig. 10. Pressure contour in Pascal in the main flow path of the valve. 




30% opening 10% opening 

Fig. 11. Contour of stream function (kg/s) inside the main valve. 



Figs. 12 and 13 depict the variation of the Mach number and the velocity vector through the valve body 
and the initial part of the spool valve. The low values of these two variables in the valve body can be attributed 
to the higher mean radius of the valve body in comparison to the bore diameter of the spool valve. As expected 
these values also remain low in the regions of the three recirculation bubbles discussed above with reference to 
Fig. 11. Between the two recirculation bubbles inside the spool valve, the flow enters the spool valve almost in 
an inward radial manner. Consistent with the decrease in the flow area due to the decrease in the mean radius, 
both the Mach number and the velocity keep on increasing. These increases signifying flow acceleration 
continue even during the flow turning to almost axial direction. Near to the axis, the acceleration persists due to 
the small inward radial component. Away from the axis up to the bubble, the flow decelerates due its diverging 
nature owing to the small outward radial component. The location of the peak velocity region away from the 
axis can be attributed to the sharper turning of the flow and the corresponding acceleration prior to reaching this 
region. Because of the complex flow structure inside the valve discussed above, CFD prediction of spool flow 
forces is more reliable. 
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30% opening 10% opening 



Fig. 12. Contour of Mach number inside the main valve. 




30% opening 10% opening 

Fig. 13. Velocity vector (m/s) inside the main valve. 

V. CONCLUSION 

Control circuit pressure characteristic has been calculated both analytically and by CFD method by 
solving compressible, turbulent flow through pressure regulated valve. In CFD analysis a special UDF has been 
developed to couple with ANSYS FLEUNT primarily to modify structured mesh of the flow domain 
dynamically with change in valve inlet pressure. CFD results has predicted better characteristic. Also, the 
simulation provides valuable insight into the control circuit behaviour, which is crucial for designers. Attempt to 
simulate using fully three-dimensional geometry is underway. 
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